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We present a method to compute off-shell effects for processes involving resonant particles at 
hadron colliders with the possibility to include realistic cuts on the decay products. The method 
is based on an effective theory approach to unstable particle production and, as an example, is 
applied to f-channel single top production at the LHC. 
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1. Introduction 

Heavy unstable particles such as the top quark, Z and W bosons as well as virtually all hypo- 
thetical particles in beyond-the-Standard-Model scenarios can only be studied through their decay 
products. The simplest way to compute cross sections involving such particles is to first con- 
sider on-shell production of the heavy particle and then its subsequent decay. Taking into account 
the matrix element of the decay of the on-shell particle — this is sometimes called the improved 
narrow-width approximation — it is possible to apply realistic cuts on the decay products. In this 
framework we can go to higher orders in perturbation theory by computing separately corrections 
to the production and the decay of the on-shell heavy particle. While this approximation is very of- 
ten good enough, this talk addresses the question on how to go beyond and include off-shell effects 
in a systematic way. Thus we are led to consider processes with resonant (nearly on-shell) rather 
than on-shell heavy particles. 

Even though the framework for computing off-shell effects presented here is quite general, 
let us consider an explicit example, ^-channel single top production. More precisely we consider 
the partonic process u{p\)b{p2) — ^ d{p^)b{p^)W^ {pw) d{p^)b{p/^)e^ {ps)v{p(,), where the 
invariant mass of the top decay products is understood to be close (but not necessarily equal) to the 
top mass, i.e. {p^ + pwY = + A with A/m^ ^ 1. The decay — )• e+ V is taken into account in 
the improved narrow-width approximation. It is clear that the dominant contribution to this process 
with its kinematic constraint is given by production and subsequent decay of a top quark. Self- 
energy resummation avoids the non-integrable singularity in the top propagator at (774 + pwY — 
by changing the denominator of the propagator to (/74 + pwY — "J? + i^t^t, where Yt is the width of 
the top. But even at tree-level, at some point we will have to take into account so called background 
diagrams which do not refer to a top quark at all. The situation becomes even more involved 
if we want to include higher-order effects. This can be done using the pole approximation 
which can be considered as a first step towards a systematic expansion of the amplitude in F^/m, 
around the complex pole. Within the pole approximation the one-loop corrections can be split in a 
gauge-invariant way into factorizable and non-factorizable corrections iQ]. Factorizable corrections 
correspond to corrections to the production and decay part of the unstable particle, whereas non- 
factorizable corrections link the two parts. Even though non-factorizable corrections (and off-shell 
effects in general) have been studied extensively in the literature they are usually neglected for 
processes at hadron colliders, because in most cases they are found to be small [^. In fact there 
are large cancellations ^ that are partly responsible for the smallness of the corrections. However, 
with cuts in the final state these cancellations are not perfect any longer and, as stated above, it is 
the purpose of this work to study these corrections. 

2. Effective theory and virtual corrections 

The salient feature in the problem at hand is the presence of different scales, rrit F,. This 
calls for using an effective theory (ET) approach to the problem. Indeed, it has been found [Q] 
that in an ET approach factorizable corrections simply correspond to the hard corrections. In 
this context, hard is understood as a mode using the method of regions [^] and means momenta 
that scale as p ~ lUt. The standard procedure within an ET approach is to first integrate out 
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Figure 1: (a) resonant tree-level diagram contributing to and A)'^^''; (b) example of a background 




(3,0). 
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diagram contributing to ^[oj^^ (c) example of a QCD diagram contributing to A'^''^'. 



the hard modes and thereby obtain an effective Lagrangian. This Lagrangian consists of gauge 
invariant operators multiplied by matching (Wilson) coefficients, which have a perturbative ex- 
pansion in the coupling and are gauge independent as well. The matching is done on-shell with 
pj = mf — imtTf. While the matching coefficients take into account the factorizable corrections, 
the non-factorizable corrections are reproduced in the effective theory by the still dynamical soft 
(and possibly other) modes. A soft mode corresponds to a momentum scaling as ~ nifS where 
we use 5 ~ A/ni,^ ~ r^/m, ~ ttew ~ (^s ^ generically denote a small quantity. The construction 
of an effective Lagrangian is a standard procedure and has been used in many different contexts. 
For the application to unstable particles it has been explained in detail in Ref. [||]. We thus restrict 
ourselves to a few remarks relevant for our example. 

The (strictly fixed-order) tree-level amplitude is a function of rrit, other masses and the mo- 
menta pi . . .p6- Choosing A = {pw + ph)^ — mf as one of the independent kinematic variables and 
expanding around the pole in A we write 
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(2.1) 



The leading part, gl^^,A^^^^y receives contributions solely from the resonant diagram, part (a) of 
Figure and scales as 5^/^, whereas for the subleading part gl^A^-^'^^ 



diagrams have to be taken into account. Given that 



- 5^/^ also background 
"^^'^ is gauge invariant it is clear that each 
term in the expansion in Eq. (2.1 ) is separately gauge invariant. Note that we also include the QCD 
contribution gewg^A^^''^^ ~ 5 which is usually considered to be a background. However, from our 
point of view this part has the same final state and contributes to £^ , albeit only at subleading order. 
Squaring the amplitude we obtain 
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This is simply an expansion of M'"'^ in the kinematic variable A which we constrain to be small. 
The (first) leading term scales as g^^j.A^^ ~ 8. The subleading electroweak and the QCD term 
scale as g^j^A^^ ~ 8^ and g^w^t ~ 8^ respectively, and omitted terms are further suppressed in 
8. This expansion is not complete yet, because it takes into account only the small parameter 
^/m^ ^ 8 but not ttgw ^Tt/nit ^ 8. The two small parameters are intrinsically linked and when 
the corresponding expansions are combined the FT becomes a very powerful tool. 
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Figure 2: Two examples of how usual loop diagrams are reproduced by the ET. Grey lines indicate soft 
gluons, solid (hollow) dots indicate tree-level (one-loop) matching coefficients. 

Within the ET framework, the leading contribution to s^^''^'^ is reproduced in three steps: pro- 
duction of an on-shell top (through an operator in the ET with its matching coefficient), propagation 
of a resonant top (through the leading bilinear top operator) and decay of an on-shell top ( through 
a decay operator in the ET with its matching coefficient). The bilinear operator resums the (gauge- 
invariant) hard part of the self energy. The subleading terms give rise to 5-point operators. Again, 
the operators as well as their matching coefficients are separately gauge invariant. 

We aim at computing the cross section up to 0'{5^l'^) ~ (^{a^d), which corresponds to in- 
cluding one-loop QCD corrections to the leading resonant part. Within the ET approach this will 
amount to computing one-loop corrections to the Wilson coefficients for the leading production and 
decay operators, as well as explicit calculation of one-loop corrections within the effective theory. 
To illustrate the connection between the ET and standard loop calculations, let us consider two 
examples, shown in Figure Q 

Example (a) illustrates how the self-energy is split into a hard part and a soft part. The hard part 
scales as 0'{8^) and, therefore, has to be resummed. The soft part is suppressed by an additional 
power 5^/^ with respect to s^^^^'^. Thus we need to include only one soft self-energy insertion 
at &{5^l'^). Example (b) illustrates how a normal loop diagram is split into a hard (factorizable) 
contribution included in the matching coefficient and a soft (non-factorizable) contribution where 
the gluon is still dynamical in the ET. The situation is not affected by including an arbitrary number 
of hard self-energy contributions to the top propagator, but an additional soft self-energy insertion 
results in a contribution beyond 

The ET allows us to identify and compute the minimal set of corrections needed at a certain or- 
der in 5. In particular a full one-loop calculation including all non-resonant diagrams is not needed 
at this order. The price to be paid for this simplification is that we have to insist on the kinematical 
constraint A/m,^ ^ 1. A corresponding calculation using e.g. the complex mass scheme would 
be more complicated, but would allow us to study observables without this constraint. 

3. Real corrections 

Using an ET approach as described in the previous section is a well established tool for the 
calculation of total cross sections [|l^]. When computing real corrections for more general observ- 
ables within an ET approach we run into difficulties. The problem is related to the fact that an ET 
approach relies on knowing and making explicit all scales of the problem. As long as we have not 
precisely defined the observable we are interested in, this is not the case. 

In Section ^ we based the expansion on the small scale A (and the couplings), implicitly 
assuming there are no further small scales. While this does put some restrictions on possible 
observables (similar to any fixed-order calculation) it allows us to study a large class of physical 
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quantities. In the case of real corrections, we have an additional gluon present in the final state. 
Depending on p-], the momentum of the gluon, the expansion parameter does change. Consider e.g. 
the case where pi is coUinear to p4, the momentum of the outgoing b quark. In this case pi and 
Pa, will combine to a jet and the proper expansion parameter would be {p/[ + pi + PwY — tnj. On 
the other hand, if p-j is well resolved, this will give rise to an additional (gluon) jet and the proper 
expansion parameter would be {p4 + pw)^ — mj, as in Section |2[ 

In order to deal with this we deviate from a strict ET approach. Using the subtraction method 
to compute the real corrections we write 

J d^n+yMn+l = J d^„+i (m„+i -M^j"^^^)) + J ^/<I>„+ 1 M^j"^^ j ^ (3.1) 

/ (m„,, -M:;- + / ^0„,iM™ 

Here M„+i is the matrix element squared for the process ub ^ dbe^ v g and denote the 

limits of the matrix elements in the singular (soft and collinear) regions. Upon phase-space inte- 
gration the term / d^n+iM^^^I^^^^ would match the infrared singularities of the full virtual one-loop 
amplitude. Since we use only the expanded virtual term, there is a mismatch between the singu- 
larities. Fortunately, the singular limits have only n-parton kinematics and we can use the same 
expansion as for the tree-level and virtual amplitudes. Computing / (i<I>„_|_iM^|"^jj'''', i.e. integrat- 
ing the expanded singular limits over the phase space, produces the same infrared singularities as 



the expanded virtual corrections. In the first term on the r.h.s. of Eq. (|3JJ) we do not perform any 
expansion. The mismatch between what we subtract and what we add back is subleading in 5 and 
beyond the accuracy we are aiming at. 



In principle we should use the full matrix element M^+i in Eq. ( 3. 1 ) to ensure gauge invariance. 
We can simplify the calculation by taking only resonant real diagrams. This potentially introduces 
a gauge dependence which, however, is beyond the accuracy of our calculation. Thus the situation 
is completely analogous to the renormalization or factorization scale dependence, which is widely 
accepted as long as it is beyond the accuracy of the calculation. A variation of /i within a reasonable 
window tells us as much (or as little) about the size of neglected higher-order terms as a variation 
of the gauge parameter § within a window around 1. Of course, it is always possible that there are 
terms that are parametrically of higher order, but numerically important. This is usually associated 
with the presence of widely different scales and requires additional resummations. 



4. Results and outlook 

We are now ready to compute an arbitrary infrared-safe quantity at As always, the 

word arbitrary has to be taken with some caution. First of all we have to enforce the constraint 
A ^ nitTt. Second, we implicitly assume there are no other small scales introduced through the 
observable. A similar requirement is present for any fixed-order calculation. 

As an example we consider ub plus jets, defining the jets using a k± cluster algorithm. 
We require a b jet, Jb, and at least one more jet, both with p± > 20 GeV. In addition we require 
150 GeV < M'''" < 200 GeV, where M'"^ = y^{pwTpjJ^ is the invariant mass of the (W Jb) pair. 
We stress that we could add cuts on the decay products of the W or any other 'reasonable' cuts. 
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Figure 3: Invariant and transverse mass distribution at ff{5) (blue, dark) and (?[5^l^) (red, light). 

In Figure | we show the M^°^ and distribution for the LHC with ^ = 10 GeV, where 
is the transverse mass of the (W Jb) system. We set = Hr = 100 GeV, m, = 171.3 GeV and 
r, = 1.32 GeV and use NLO MSTW pdfs ^ for the LO and NLO results. Compared to previous 



NLO calculations 012| ] our result also includes non-factorizable corrections which leads to a small 
but visible deviation from a Breit-Wigner shape in the in the M'"^ distribution. 

Of course, the results shown in Figure ^ are not complete, since at NLO we also have to 
include gluon-initiated processes and there are other partonic initial states such as bu ^ dbW and 
cb ^ sbW . A full analysis as well as the potential impact of non-factorizable corrections on a 
measurement of nit for single top and in particular tt pair production is left for future work. 
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